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ABSTRACT 



We solve the one loop effective scalar field equations for spatial plane waves 
in massless, minimally coupled scalar quantum electrodynamics on a locally 
de Sitter background. The computation is done in two different gauges: a 
non-de Sitter invariant analogue of Feynman gauge, and in the de Sitter 
invariant, Lorentz gauge. In each case our result is that the finite part of the 
conformal counterterm can be chosen so that the mode functions experience 
no significant one loop corrections at late times. This is in perfect agreement 
with a recent, all orders stochastic prediction. 
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1 Introduction 



Gravitons and massless, minimally coupled (MMC) scalars are unique in 
possessing zero mass without classical conformal invariance. This allows 
them to mediate vastly enhanced quantum effects during inffation [1]. Many 
examples have been studied over the history of inflation, starting with the 
first work on scalar [2] and tensor [3] density perturbations. Those were tree 
order analyses. More recently, there have been analyses of loop corrections to 
density perturbations [H El [6l [7] and to similar fixed-momentum correlators 

A reasonable paradigm for much of inflation is de Sitter background and 
a wide variety of enhanced quantum loop effects have been studied on this 
background in many different theories. In pure quantum gravity the one loop 
self-energy [9j and the two loop expectation value of the metric [TOl [11] have 
been calculated. For a MMC scalar with a quartic self interaction the expec- 
tation value of the stress tensor [121 US] and the self- mass-squared [H] have 
both been computed at two loop order. In scalar quantum electrodynamics 
(SQED) the vacuum polarization [151 [IB] and the scalar self-mass-squared 
[TTj were evaluated at one loop order, and the expectation value of two scalar 
bilinears have been obtained at two loop order p3] . In Yukawa theory the one 
loop fermion self-energy [TSl [20] and the scalar self-mass-squared [21] have 
both been evaluated at one loop order. A leading logarithm resummation 
has also been obtained for Yukawa, and checked against an explicit two loop 
computation [22]. And the one loop fermion self-energy has been calculated 
in Dirac + Einstein |23j . 

Of course density perturbations are directly observable, as are the expec- 
tation values of operators such as the stress tensor. However, an additional 
step is necessary in order to infer physics from a one particle irreducible (IPI) 
function such as the vacuum polarization [^11''] (x; x'), the fermion self-energy 

These IPI functions 



[iTjj]{x; x'), or the scalar self- mass-squared M {x;x' 
correct the classical linearized field equations as follows. 



{x;x')A^{x') = 



d,{V^g''Pg^-J'p^{x)) + J d'x' 

^^[^/^g^''^^<l>{x)) - J d^x'M\x;x')<!>{x') = . 
One can infer physics from the quantum-corrected mode functions. 



(1) 
(2) 
(3) 
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Doing this for photons in SQED [TT], and for fermions in Yukawa [H], 
leads to the surprising result that the inflationary production of scalars en- 
dows these particles with mass. For fermions in Dirac + Einstein one finds 
that the field strength grows by an amount which eventually becomes non- 
perturbatively large [21]. On the other hand, the result for the scalars of 
Yukawa theory is that the finite part of the conformal counterterm can be 
chosen so that there are no significant late time corrections at one loop order 
[21j. The purpose of this paper is to demonstrate that the same result applies 
for the scalar mode functions of SQED. 

We begin in section 2 by defining the effective mode equation and the 
sense in which we solve it. We also explain the relation between the €-number 
solutions we shall find and the quantum operators of SQED. Of course the 
self-mass-squared of SQED is gauge dependent, and has been computed in 
two different gauges. Section 3 solves for the mode functions in a non-de 
Sitter invariant version of Feynman gauge [251 [U]; section 4 does the same 
for the de Sitter invariant Lorentz gauge jMlETj. In both cases our result is 
that the finite part of the conformal counterterm can be chosen so that there 
are no significant one loop corrections at late times. Section 5 summarizes 
our work and discusses its implications. 

2 The Effective Mode Equation 

It turns out that the same operator formalism gives rise to many different 
effective field equations in quantum field theory. The purpose of this section 
is to specify both the particular ones we are solving and the sense in which 
we shall solve them. We begin by recalling the operator formalism of SQED 
[T71 ri8j . We then parameterize general quantum corrected mode functions 
and contrast the parameter selections appropriate to fiat space scattering 
problems with the choices appropriate to cosmology. A brief review is given 
of the Schwinger-Keldysh formalism whose linearized effective field equations 
give the quantum corrected mode functions for cosmology. We close with a 
discussion of the limitations on our knowledge and what they imply about 
the sense in which we should solve the effective mode equation. 
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2.1 Relation to Fundamental Operators 



We work on the nondynamical metric background of de Sitter in conformal 
coordinates with spacehke signature, 



Quudx^dx" — a^(r))(—drj^ + dx-dx) where a(rj) — ——— . (4) 
^ ^ Hr] 



In other words, (jf^j, = a^rj^y, where 7]^^ is the Minkowski metric. When 
expressed in terms of renormahzed fields and couphngs, the Lagrangian of 
SQED takes the form. 



-SZ2{d, 



V9 



--5Z^F^F^,gP^g° 



-9 - S^^*^R\^ 

r^-hx{ip*ipr^. (5) 



We have chosen to study the exactly massless version of the theory in which 
the renormahzed values of the conformal and quartic couplings vanish. Hence 
the parameters have the following expansions in terms of the charge e. 



= e^6Z20 + 6^5^2,4 



5X = e^5Xi + ... , 
5Z^ = e^6Z.2 + e'^SZ 



3,4 



(6) 
(7) 



The corrected mode functions for which we will solve are related to the 
fundamental operators through a number of choices which require extensive 
explanations. To fix notation for these explanations we begin by stating the 
relations. 



k) = 
A^i{x; k, A) = 



ip{x),a\k)]\'^i 
A^{x),^\k,X) 



(8) 
(9) 



The quantities to be explained are the states and |^'/), and the free 
creation operators a'^^k) and 7^(A^, A). 

Flat space scattering problems correspond to taking to be the state 
whose wave functional is free vacuum at asymptotically early times. One 
also chooses to be the state whose wave functional is free vacuum at 
asymptotically late times. Although these choices have great physical interest 
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for flat space scattering problems, they have little relevance for cosmology. 
In cosmology the universe often begins at a finite time, and it evolves to 
some unknown state in the asymptotic future. Persisting with the choices of 
flat space scattering theory would result in acausal effective field equations 
which are dominated by the imperative of forcing the universe to approach 
free vacuum. The fact that we are looking at matrix elements, rather than ex- 
pectation values, would also have the curious consequence of giving complex 
results for the matrix elements of Hermitian operators. 

In cosmology it is better to imagine releasing the universe from a pre- 
pared state at some finite time, and then watching it evolve as it will. This 
corresponds to the choice |\E'/) = |\E'j). We additionally assume that both 
are free vacuum at ri = rii. 

So much for the states in relations ([8ll9]). To understand the free creation 
and annihilation operators one integrates the invariant field equations of 
SQED, 



1+5 Z; 



+5Z2 2(1+5^2) 



0, 

: . 



(10) 
(11) 



Here the covariant derivative operator is = d^ + ieA^. The result of 
integrating (llOtllip is the Yang-Feldman equations [28] , 



A^{x) = Ao^{x) + 
The two interactions are, 



Vi 



dri I d^~^x'G{x;x')I[ip\ip,A]{x') 



'12) 



d^'-'x' 



{x;x')r[ip\ip,A]{x') . (13) 



I[ip*,ip,A] 



jnV*,^,A] 



-ieA,^ 



gg^'d,V-ied,{^gg^''A,^) 

5i\/—gR^ (5 A A 



—le 



+V^g'"'e^A^A,if 

'gg'"' 



-gif*ip^ 



1+5Z2 



2(1+5^2) 



1 + 5^3 

G{x\ x') is any solution to the free scalar Green's function equation, 

d^[-^g^"'duG{x;x')) = 6^{x-x') , 



(14) 
(15) 

(16) 
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and [^iGi] (x; x') is any solution to the analogous photon Green's function 
equation in whatever gauge is employed. For flat space scattering problems 
the Green's functions would obey Feynman boundary conditions — and hence 
amount to —i times the Feynman propagators. In cosmology it is more 
natural to use retarded boundary conditions. 

We stress that the fundamental operators '{){x) and A^{x) are unique and 
unaffected by the choices of r]i and the boundary conditions for the Greens 
functions. What changes as we vary rji and the Green's functions is the free 
fields, ^q{x) and y4o^(x). The fact that they obey the linearized equations 
of motion and agree with the full fields at t] = •)]{ implies that they can be 
expanded in terms of free creation and annihilation operators, 

Mx)= /(^^{^(^' k)e'' Mk) + ^*iv, k)e-''-'P\k)} , (17) 
Ao,{x)= k, A)e^S(fc, A) + e;(r/, k, X)e-'^'-'j\k, A)} . (18) 

Here u{ri, k) and 6^(77, k, A) are the Bunch-Davies mode functions [29j in what- 
ever gauge is being used. They also have canonically normalized Wronskians, 
which for the scalar is, 

u{r], k)doU*{r], k) - dou{r], k)u*{r], k) = —— . (19) 

Although the various creation and annihilation operators change as different 
Green's functions are employed in (fT2!) and (fT3|) . their nonzero commutation 
relations remain fixed, 

[«(fc),«t(fc')] = (27r)^-i(5^-i(fc-^') = [/?(^),/3^(^')] , (20) 
7(fc,A),7^(^',A')] = (27r)^-i5^-i(^-^')5AA' . (21) 

One develops perturbation theory by iterating (11211131) to obtain expan- 
sions for the full fields, '^[x) and A^{x\ in terms of the free fields, ^q{x) and 
Ao^(x). With our choice of retarded boundary conditions the expansion for 
'^{x) takes the form, 

ip{x) = Mx) + f'dr]' [ d''~'x'G,,tix;x')I[ip\ip,A]{x') , (22) 

= Mx) + r dv' I d^'-'x'Gretix; x')I[ipl, ipo, Ao]ix') + ... . (23) 
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With our choice of = l^^i) as free vacuum at rji, the quantum corrected 
scalar mode function is, 



A;) 



u(r7,A;)e*'=-^ + 0(e^) 



(24) 



2.2 The Schwinger-Keldysh Formalism 

It is straightforward to demonstrate [24j that the quantum corrected mode 
function 0241) obeys the hnearized effective field equations of the Schwinger- 
Keldysh formalism [501 [5T| [23] • This is a covariant extension of Feynman 
diagrams which produces true expectation values rather than the in-out ma- 
trix elements of conventional Feynman diagrams. Because many excellent 
reviews of this subject exist [311 ESI EHl EH H] we will simply state that part 
of the formalism which is necessary for this work. 

The basic idea is that the endpoints of propagators acquire a ± polar- 
ity, so every propagator iA{x; x') of the in-out formalism generalizes to four 
Schwinger-Keldysh propagators: iA+_|_(x; x'), iA+_{x; x'), iA_+(x;x') and 

iA {x;x'). The usual vertices attach to + endpoints, whereas — endpoints 

attach to the conjugate vertices. Had the state been some perturbative cor- 
rection of free vacuum, this correction would give rise to additional vertices 

on the initial value surface h ones from 'ifi[ip^{rii),ip^{rii)] and — ones 

from 'iff[ip*_{r]i),(p_{r]i)]. 

Because each external line can be either -|- or — in the Schwinger-Keldysh 
formalism, each A^-point IPI function of the in-out formalism corresponds to 
2^ Schwinger-Keldysh A^-point functions. The Schwinger-Keldysh effective 
action is the generating functional of these IPI functions, so it depends upon 
+ fields and — fields. The effective field equations come from varying with 
respect to either polarity and then setting the two polarities equal. 



6(pl{x) 



d7]' j rfV{M^^(x; x') + MlSx] x')]^{x') + 0(y?>2) . (25) 



Note that we have taken the regularization parameter D to its unregulated 
value of D = 4, in view of the fact that the self-mass-squared is assumed to be 
fully renormalized. It is the linearized effective field equation which $(x; fc) 
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obeys, 

(26) 

Converting the single, in-out self-mass-squared — M^(x;x') — into the 
four polarities of the Schwinger-Keldysh formalism is simple. The +-|- po- 
larization is the usual in-out self-mass-squared and the polarization is 

minus its conjugate. 



Ml^{x;x') = M^{x;x') 



and 



ML (x; x') = - [m^ (x; x')) * • (27) 



To obtain the mixed polarizations we observe that, in de Sitter conformal 
coordinates, the in-out result depends mostly upon the conformal coordinate 
interval |24j , 



Ax^_^(x; x) 



x—x 



ri — ri'\—i5] . 



(28) 



To obtain the H — and 
things 



polarities at one loop order we must do three 



Replace AxL(x;x') with the appropriate coordinate interval, 

2 / , \ 2 



Ax^_(x;x') = 
Ax^^(x;x') = 



x—x 



— [rj — ri' + id 

— {rj — rj' — i5^ 



(29) 
(30) 



• Drop all delta function terms; and 

• Multiply the result by —1 to account for the — vertex. 

It will be seen that the ++ and H — terms in fl26l) exactly cancel for rj' > r) 
and also, in the limit 6^0, for x'^ outside the light-cone of x'^. This is how 
the Schwinger-Keldysh formalism gives causal effective field equations. 



2.3 The Meaning of "Solve" 

Two important limitations on our current knowledge restrict the sense in 
which we can usefully solve the effective mode equation. The first of these is 
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that we do not know the exact scalar self-mass-squared. The full result can 
be expressed as a series in powers of the loop counting parameter e^, 



M++{x;x') + M+_{x;x') = Y^e^^ M\[x]x') . (31) 



£=1 



Because we possess only the l= \ term, we can only solve the effective mode 
equation (126!) perturbatively to order e^. That is, we substitute a series 
solution, 

oo 

^{x- k) = uiri, k)e'^'^ + J2 e^^^iiv, k)e'^'" , (32) 

1=1 

and then segregate according to powers of e^. The 0th order solution is well 
known, 

/ 7 \ r 1 r 1 / \ 

^W'^) = —fKf^ 1~T}- exp — ,. (33) 



'2k^ 

The 1st order solution $1(77, k) obeys the equation 



dl + 2Hado + k'^]^i{r],k) = - / di]' j d^x' Ml{x]x')u{ri\ky^<^-'^"^ . 

(34) 

The second limitation is indicated by the subscript on the temporal 
integration in we release the universe in free vacuum at time ri = r]i. Not 
much is known about the wave functionals of interacting quantum field the- 
ories in curved space, but it can hardly be that free vacuum is very realistic. 
All the finite energy states of interacting flat space quantum field theories 
possess important corrections. It is inconceivable that similar corrections are 
not present in curved space, at least in the far ultraviolet regime for which 
the geometry is effectively flat. 

Although one could perturbatively correct the free states to a few orders, 
just as in nonrelativistic quantum mechanics, the standard procedure of fiat 
space quantum field theory is to instead release the system in free vacuum 
at asymptotically early times. In the weak operator sense, infinite time 
evolution resolves the difference between free vacuum and true vacuum into 
shifts of the mass, field strength and background field [38j. In cosmology 
we cannot typically employ this procedure, however, it is still possible to 
perturbatively correct the state wave functionals. 

Corrections to the initial state would show up as new interaction vertices 
on the initial value surface. One would expect them to have a large effect 
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on the expectation values of operators near the initial value. One would also 
expect their effects to decay in the expectation values of late time operators. 
We have not worked out these surface vertices but the need for them has 
been apparent from the very first fully regulated computations of this type 
[T2I [T6] . For example, consider the expectation value of the stress energy 
tensor of a massless, minimally coupled scalar with a quartic self interaction 
on a locally de Sitter background. If the state is released in free vacuum at 
the instant when the de Sitter scale factor takes the value a = 1, then one 
can choose the renormalization parameters so that the energy density and 
pressure induced at two loop order are [T3] . 

P = 




P = 

It is the terms which fall off like powers of 1/a that we suspect can be 
absorbed into an order A correction of the initial state, 

= i^A~lt^-^^^\+0(X^), (37) 





AP = 7^^^-i:^^^y^K0(A^). (38) 

Of course the fact that they fall off as one evolves away from the initial value 
surface suggests that they can be absorbed into some kind of local interaction 
there. Note also that they are separately conserved, which is exactly what 
would be the case if they could be canceled by a new interaction vertex. 

Note particularly that expressions fl37ti38p diverge on the initial value 
surface. Divergences on the initial value surface have also be found in the 
effective mode equations for photons in SQED [I6l[39|, for fermions in Yukawa 
theory [19], for scalars in Yukawa theory [21], and for fermions in Dirac 
+ Einstein pi]. These initial value divergences reflect the fact that free 
vacuum is very far away from any physically accessible state. Equation 
accurately determines the one loop correction to the mode function ([24 
appropriate to free vacuum, however, that mode function has little physical 
relevance because free vacuum can not be assembled. 

If one desires expressions for physically relevant mode functions which are 
valid even for times near the initial value surface, there is no alternative to 



9 



including corrections to the state wave functional in the self-mass- squared. 
This seems quite practicable because one need only do it perturbatively to 
the same order as the ordinary, "volume" contributions are known. Although 
there are no stationary states for SQED on de Sitter background, it would 
probably be an excellent approximation — becoming exact in the ultraviolet 
— to simply solve for the relevant corrections in flat space. 

Unfortunately, we do not now possess the order A correction to the state 
wave functional. That means there is absolutely no point in trying to solve 
for all times. However, it does not mean (l34ll is devoid of physical 
information. In particular, the effects of the state corrections must fall off 
at late times — and quite rapidly, like inverse powers of the scale factor. 
This fall-off is evident in expressions ( 1571155]) . It simply reflects the same 
process which is employed in flat space quantum field theory [38] whereby 
time evolution washes away the difference between free vacuum and true 
vacuum. We can therefore extract valid information from (|34|) by solving it 
in the late time limit. 



3 Non-de Sitter Invariant Gauge 



Reliable results for the one loop scalar self-mass-squared have been obtained 
in two gauges. The first is a non-de Sitter invariant analogue of Feynman 
gauge which corresponds to adding the gauge- fixing term [17], 

1 



GF 



,D-4 



(39) 



In this gauge the renormalized in-out (and hence also, the -|-+) self-mass- 
squared is. 



-g6'^{x — x') 



e^aa' 



871^ 



ln{aa')d'6\x-x') 



■ \n{a)5'^{x — x') 



ie^aa 



'-d4\n^(^H^/\xl^)-2\Ti(^H^/\x] 



[aa 



/\3 



In^ 



V4 +- 



92 r 



ln^(Ji7^Axg-31n(^i7^Aa:g] + O(e^) . (40) 
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Here e^SZ^"^ and 6^5^^°"^ are the arbitrary finite parts of SZ2 and S^, in tliis 
gauge and at order e^. Tlie various derivative operators are, 

n = -^d^{y^g''''d,) , d''=r]^''d^d, and = dA ■ (41) 

Applying tlie rules of the previous section allows us to recognize the + — 
self-mass-squared as, 



Mnon {X; X') 



2%4 I V4 +-J V4 ^ 



\n^(-H^Axl_)-ln(-H^Ax^ 



'\n''{^H^Axl^) -3\n{^H^Axl_)] | + O(e^) . (42) 

The next step is to combine fHOj) with P2|) to read off Ail{x] x') for the 
right hand side of fl34l) . To simplify the notation we first define the spatial 
and temporal intervals, 

Ar = — x'll and Arj = rj — rj' . (43) 

We also take note of the differences of powers of ++ and H — logarithms, 

\Yi(^H^Axl^)-hi(^H'^Axl_) = 2me{Ar]-Ar) , (44) 

\n\^H^Axl^)-\n\^H^Axl_) = Ame{A7]-Ar)\n[^{Ar]^-Ar'')]. (45) 
With this notation we find, 

— - \Q.{a)5^{x~x') 

471-2 \ J \ I 



+ ^ ln(aa')a25 



ad 



267r3 



- d^\e{Ari~Ar) {^^{Aif - Ar'^~^-\\^ | 
+ (aa')4 V ^(Ar^-Ar) f Inf— (Ar^'-Ar') 



r 



92 



/ tj2 

^ ( A77 - Ar ) I In [— ( Ar]^ - Ar 2 ^ 



3^ 



>. (46) 
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Note that the result is manifestly real and that each of the nonlocal factors is 
shielded by a causality-preserving factor of ^(At^— Ar). These are features of 
the Schwinger-Keldysh formalism which are absent from the usual (in-out) 
effective field equations. 

We now evaluate the contribution each term in fH6l) makes to the right 
hand side of the one loop effective mode equation flMl) . The local contribu- 
tions are simple, 



rO 



u{ri, k)e'^-'^ = , (47) 



ik-x 



-[-AMHiZ'u{r,,k)} , (48) 



\'[d^x'\ ^ \n{aa')d^6\x-x') u{r]', k)e'^-^'^'-'^^ 

,; J OTT 



^l\n{a){dl + e){au{r],k)) + [d^, + e)[a\n{a)u{r],k))\ , (49) 



Stt 

H^a^ f . , . . . , . 3 , , , 1 



2 ln(a)M(?7, A;) + -^(r], k) + —dou{ri, k) ) , (50) 



47r2 [ ' ' 2 ' Ha 
- [\'[d^x'l-^aHn{a)6\x-x') u{r]', k)e"^<'^'-'^^ 



Vi 



= ^{Ha)uiv,k)}. (51) 

The contribution of the first nonlocal term in fH6l) has been evaluated 
in a previous study of one loop corrections to the scalar mode functions of 
Yukawa theory [21j. By making the replacements, 

f — ^ 1 , ^ — ^ y and g{r]', k) — ^ u{t]\ k) , (52) 

we can read off the result from equation (68) of that paper, and some of the 
subsequent asymptotic expansions, 

-e-*"^aV Va'M(r/', A;) / e'^""|lnf-i/2(Ar/2- Ar^) 

' Jn JAr<Av '-4 



2^77^ Jrii JAr<Ari 
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a 



47r 



{—a \ia(a)u(r], k) 
+do i:;jv'a'uiri',k)cos{kAr])\nil-^) \ ,(53) 
+kj^^ dri'a'u{T]', k) sm{kAri) ln(l-f ) 

2\n{a)u{r],k)-3u{r],k) + 0[\) \. 



(54) 



We evaluate the contribution of the second nonlocal term in (H6|) by first 
making the change of variable f=x' — x, then performing the angular inte- 
grations, and finally changing the radial variable to r = Ar]z, 



4^3 



47]-2 



^~^^^■xyl dr]'a'\{r]',k) / d'x' e 

JArKArj 



ik-x' 



In 



rl 



H'^{Arf-Ar'^ 



k jdT]'a'^u{T]\ k)ATfjdz z sm{kAr]z)l2 ln(i7 Ar^) + ln(i^) - ^| .(55) 



The z integration can be performed in terms of sine and cosine integrals but 
there is no need to do this. We simply observe the behavior of the z integral 
for small A?7, 

j dzz sin(A;Ar/z) |2 ln(if Ar]) +ln(^^) - ^ | — > '^kArj \n{HAr]) . (56) 

This implies that the rj' integration in fl56l) converges at arbitrarily late times 
{i]—>-0^). Hence the second nonlocal term can be expanded as follows, 



4^3 



I£a 



e-'^-^V^rdr]'a"^u{7]\ k) / d^x' e 



Vi 



j3^' Ak-x' 



In 



Ar<Ar] 



H^{Arf-Ar' 



4 

47r2 



A- + o(!^)[, (57) 



where the constant K is, 
k r°° /ui^r]', k) 



K 



n^'^^^ f'dz z sm(^) I -2 In(a') + lnf 
Jai a' Jo ^Ha'^ ^ 



l-z\ V 



(58) 



We will see that only terms which grow as fast as can give significant 
corrections at late times. 

The contribution from the final nonlocal term in fH6l) is the most difficult 
to evaluate. The initial steps are the same as for the second term. However, 
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one must then act the temporal derivatives and express the z integration in 
terms of sine and cosine integrals, 



4^3 



2%^ 



-'^■^d'^ I "dr]'a'\{7]', k) I d^x' e'^-"^' < 



In 



rl 



Ar<Ari 



4^3 



8^k 



{d'o+k') dr]'a"u{7]',k)A7] 



X dz zsin{kAr]z)<2\n{HAr]) + \n 



I: 



Sir^k 



2a 



1-z' 



3 



3] 
21 



(59) 



■ [dri'a'^u{r]' , k)\ —2 cos{a) [ dt 
Jrii [ Jo 

+2 sin(a) 



sin(t) 



2°, cos(t)-l 
dt — 



-21n(^)-l 



(60) 



Here we define a = kArj. 

To extract the leading late time behavior of expression fl60l) we first note 
that the term in curly brackets vanishes like 4aln(Q;) for ri' = r]. This means 
that the upper limit makes no contribution if we integrate by parts on a'^, 



'dv'a'Mv',k){ } = ^u{v',k){ } 



1 n 



2H Jr,^ 



di^'a'^d',[u{r^\k)[ }],(61) 



■^a]u,{ }^-^J\'a"d',u{v',k){ }+u{v',k)d',{ } 



■ (62) 



Of course the term from the lower limit approaches a constant at late times, 
and it need not concern us further. The same is true of the integral in fl62!) 
which involves the derivative of «(//', k), 



dou{r]', k) 



H 



Ha' 



r ik 



exp 



Ha' 



(63) 



That leaves only the integral which involves the derivative of the curly brack- 
eted term of fl60l). 



2° cos(t)-l ,,,,,, 1 
dt +2 ln(i7Ar/) — 



&q[ } = -2k sin(a) j^'^dt ^^^-2k cos(a) 



(64) 
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The term we need comes from substituting fl64l) into the final integral 
of expression fl62l) . The contribution from the two t integrals approaches a 
constant at late times, 



k /■'',/ ,9 //, X /■ , rCOs(t-Q;)-cos(a) 
— dr]' a'^u{7]\k) / dt ^ ^ ^ ^ 

H Jrii Jo L 



t 



constant . 



(65) 



Each of the two remaining terms in flMl) requires a further partial integration 
to extract the leading late time behavior. For the term k cos(a) we integrate 
by parts on a'^, 

k /"'? 

— : / di]' a''^u{ri' , k) cos{a) 

2H Jr]i 



2m 

k 
2IP 



auiji' , k) cos(a) 



au{ri, k) + constant . 



' r ~\ 

drj' a' O'q u{ri' , k) cos{a) , (66) 

(67) 



The logarithm requires that we partially integrate on a''^hi{HArj) using the 
relation, 



fdr]'a'^ln{HA7]) = -^\{a~a') ln(l--)+a'ln(a') 



(68) 



The result is, 
2k n. 



2k 2k 
— 77 / drj' a'"^ \n{H Ar])u{r]' , k) cos(q;) = —77^ ahi[a)u{7], k) + constant . 

(69) 

It follows that the late time expansion of the final nonlocal term in (H6l) has 
the form, 



4^3 



ma 



-ik-x p\2 



d'l d7]'a'%{v',k) / d^x'e 



3^/ Ak-x' 



In 



-^H\Ar]^-Ar^] 



^^-^|- ln(a)M(r/, k) + ^u{r], ^) + 0(^ 



2] 
(70) 



Summing the various local and nonlocal contributions, and recalling that 
dou{r],k) ~ —k'^/Ha x u{0,k), we find that the one loop effective mode 
equation in this gauge takes the form. 



dl+2Had^k^ ^f^'l^, k) 



2^4 



ma 

47r2 



^87r25ar]«(^,^)+0(i)|. (71) 
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To infer the asymptotic solution it is best to change the derivatives from 
conformal time t] to co-moving time t = — \n{—Hr]) / H , 



dl + 2HadQ + A;^ 



(72) 



By choosing the finite part of the conformal counterterm appropriately, 

5 1 



non 
fin 



12 (47r) 



(73) 



we reach an equation of the form, 



(74) 



where C is a constant. Except for possible homogeneous terms — which can 
be absorbed into the finite part of the field strength renormalization at late 
times — the one loop correction rapidly redshifts to zero. 



C 1 



(75) 



4 de Sitter-Lorentz Gauge 

Reliable results also exist for de Sitter-Lorentz gauge [26| 127]. 



(76) 



Because this gauge preserves manifest de sitter invariance it is best to express 
results in terms of the following function of the invariant length i{x; x'), 



y[x;x 



4sin^ 



-Hi{x;x') =aa'H\x-x'Y{x-x'yr]^^ . (77) 



The imaginary parts appropriate for the polarizations of the Schwinger- 
Keldysh formalism are, 

y++{x;x') = aa'H'^Ax'^_^_{x]x') and y+.{x;x') = aa' H'^Ax'^_{x; x') . 

(78) 
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In this notation the scalar self- mass-squared is 



M inv I ^ . 



-e bZ^^ . 



(An) 



-g'n 



+ 



(47r)4 ■ 
6 

1 y++ 
^ 4 



4 ' 



Before giving the H — term and combining to work out Ai1{x;x'), it is 
best to remove the factors of 1/y using the identities, 



4 

y 



Ki) - ^{i>"Hi)-m))+3m)-2. (81) 



Because any analytic function of y{x; x') will cancel in the sum of ++ and 
H — polarizations, we do not write out such terms explicitly. 



+ 



+ 



(47r)4 
□ 



V 4 



81n(l 



^ln2f^)+ln( 



+ 



i74 



21n2f^ 
V 4 



1/4 



+AnalyticV + 0(e*) . (82) 



The procedure of Section 2 gives the corresponding H — self-mass-squared. 



M?nv (x; x') 



(An) 



-9\/-9'< 



-iln^fcUlnfc 



+ ^4 



2.„=(^) 



□ 



-41n^(^)+81n(l-^)ln(^) 



+121n(^)-^"^'^ 
V 4 ^ 1 



4 



+ Analytic 



4 



(83) 
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Of course the delta functions terms in Ai1{x;x') make the same contri- 
butions as they did for the noninvariant gauge of the previous section, 

- j\'Jd'^x'{-5Z'^Ja''n5\x-x') + l2H''5^Za^S\x-x')}uir]\k)e^^^^^^ 

= -12H^6^fina\{v,k) . (84) 

We can avoid a lengthy computation of the contributions from the nonlocal 
terms by taking note of four points: 

• As explained at the end of Section 2, the only sensible and physically 
interesting regime in which we can solve the effective mode equation is 
for late times, long after the state was released; 

• In this regime we can replace the 0th order mode function by a con- 
st ant 

- rdri'fd'^x'Mlix; x')u{7]', fc)e*(""-"^) ^ -u(0, k) (\'(d'^x'M\{x- x') ; 

(85) 

• This reduces the nonlocal contributions to a sum of terms of the form, 

• The four integrals of the form ( l86l) which we require have been worked 
out in a previous computation of the vacuum expectation value of two 
coincident scalar bilinears at two loop order |18j. 

Table [1] reproduces the necessary integrals from previous work ^18] . Be- 
cause the integrals can only depend upon a, the following identity is useful 
for acting the d'Alembertians, 

(a ln2(a) + (3 ln(a) + 7 + 5^^) = -a f 6 ln(a) + 2) - 3/5 + 5 T^M^" . 

(87) 



^For u{ri' ,k) this is obvious from the fact that expression (j33[) approaches a nonzero 
constant at late times. That the spatial plane wave factor of e*'^ ^^ can also be dropped 
follows from the causality of the Schwinger-Keldysh formalism. The factors of ^(ATy— Ar) 
which arise whenever H — h and H — terms are added — for example, expressions (j44ll45p — 
require — a;|| <i] — ri', which is small in the late time regime. 
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fix) 


jdVa'H/(^)-/(^)} 


\n(x) 
1-x 




ln2(a;)-21n(l-a;) ln{x) 


1 7r2 1 ^/ln(o)\ 
6 9 ^ a ^ 


ln{x) 


iln(a)-i+0(i) 




|ln^(a) - |.n(a) + J - f + OC^) 



Table 1: Integrals of de Sitter Invariants. Multiply each term by ^^t^} . 



The triple d'Alembertian nonlocal contribution is, 

^3 ^0 r 1 



- „(o / -1^ + 1 '"(«) - i + S 1 (88) 

-"(°'*' 16.^ H'\ +ll„(a)-ll + 0(!=M) }• 

= »(0,*)^|0 + O(^)|. (89) 

The double d'Alembertian term gives, 

^ ' ^ 167r2 /f4|3 V ; g V y 2 9 V a ^J' ^ ^ 
= ''(0.'=)15^{6 + O(^)|. (91) 

And the term with only a single d'Alembertian is, 

(n t^!^^!^-^/?/ ' '4/^3 '/-41n'(^)+81n(l-^)ln(^)-(+-)1 
""^"'^^ (47r)4 W.r'' y''^l+121n(^)-61n(^)/(l-^)-(+-)/ 

"''^^'''^ 167r2 +21n(a)-^-3 /' ^^^^ 

,,3i7V r .ln(a)xl 
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Summing the local contribution fl84l) . and the three nonlocal ones fl89p . 
(^T]) and (19^ . gives the following effective mode equation in Lorentz gauge, 



d^ + 2Hado + k^ 



2^4 



167r2 



127r25eiiXr7,A;) + 



ln(a) 



(94) 



This is quite similar to the result of the previous section, and the subsequent 
analysis is the same. By choosing 5.^^°'' = 0, and converting to co-miving 
time, we obtain the form. 



d^ + 3Hdt + ^Wriv,k) 



167r2 



C 



ln(a) 



(96) 



where C is a constant. Except for possible homogeneous terms — which can 
be absorbed into the finite part of the field strength renormalization at late 
times — the one loop correction rapidly redshifts to zero. 



C ln(a) 



327r2 



+ 0(1). 



(96) 



5 Discussion 



The Schwinger-Keldysh formalism gives effective field equations which are 
causal, and which allow real solutions for Hermitian fields. They are well 
adapted to cosmological settings in which the physically sensible experiment 
is to release the universe in a prepared state at some finite time. If this 
state is chosen to be free vacuum, the formalism will be simple in the sense 
of lacking interaction vertices on the initial value surface. However, the 
unphysical initial state results in effective field equations which diverge on 
the initial value surface. Even at late times, there will be some contamination 
from the unphysical initial state in the form of terms which decay. This could 
be avoided by perturbatively correcting the initial state — a worthy project 
which we have, unfortunately, not undertaken. However, one can still employ 
the simple equations reliably at late times where there are no divergences and 
the contamination is dying away. 

We have solved the free-vacuum equations for one loop corrections to 
the scalar mode functions of SQED on a locally de Sitter background. Be- 
cause the mode functions are gauge dependent, the computation was made 
in two different gauges. In each case it was possible to choose the finite part 
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of the conformal counterterm so as to prevent the occurrence of significant 
corrections at late times. 

This might seem a very different outcome from similar studies of one loop 
corrections to the mode functions of photons in SQED [16l [39], fermions in 
Yukawa theory [191 120] and fermions in Dirac + Einstein [23]. In each of 
those cases the mode functions suffer one loop corrections which grow at late 
times, instead of falling off. In fact our result for the scalar mode functions of 
SQED does fit this pattern for a generic value of the conformal counterterm. 
Had we not chosen 6^^^ = in (IMI) the one loop correction would obey the 
equation, 

'd^+3Hdt+^]^, = -lH'5(,M0,k) + O(^-^) . (97) 
In that case the resulting solution would indeed grow at late times, 

$1 = -^56nu(0, k) ln(a) + O(i^) . (98) 

So the difference between the scalar of SQED and the other cases is just that 
SQED possesses a free parameter which can be tuned to suppress the large 
one loop corrections that would otherwise occur. 

It is worth pursuing this point a little further. Based upon fl98l) . the form 
one expects for the largest late time corrections to the full mode function is, 

oo M ( °° £^ 

<!>{^,k) =uiv,k) + Y.e''Mv,k) ^\l + ^ce{eHn{a)) . (99) 
1=1 v^fc [ J 

Secular corrections which involve powers of the infrared logarithm ln(a) are 
ubiquitous in quantum field theories that include MMC scalars and/or gravi- 
tons [T0l[IIl[ll[T3ll31[231[8l[Tg. The continued growth of ln(a) offers the 
fascinating prospect of compensating for the small loop counting parameters 
— in this case — which usually suppress quantum loop corrections. How- 
ever, the valid conclusion from a series of the form (l99l) is that theory breaks 
down at ln(a) ~ 1/e^, not that quantum loop effects necessarily become 
strong. A nonperturbative analysis is required to determine what actually 
transpires. 

Starobinskii has developed a stochastic formalism [IQ] which has been 
proven to reproduce the leading infrared logarithms at all orders in scalar 
potential models [H] |12] . Starobinskii and Yokoyama have shown how this 
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formalism can be used to obtain explicit, nonperturbative results for general 
expectation values in these models |13]. Starobinskii's formalism has recently 
been extended to Yukawa theory [23] and, of special interest to us, to SQED 
[Hi |l5]. This nonperturbative formalism allows us to view the results we 
have obtained in a larger context. 

First, note that the choice of 5.^^°^ = 0, which enforces the absence of sig- 
nificant late time corrections for de Sitter-Lorentz gauge, coincides precisely 
with the choice which cancels the leading infrared logarithm corrections to 
the expectation value of ip*{x)ip{x) at order in the same gauge [18]. This 
same choice cancels the (p*ip contribution to the effective potential [HI B5] . 
The concurrence of three distinct results seems to confirm the stochastic pre- 
diction [m US] that significant late time effects on the scalar derive solely 
from the effective potential and not, for example, from corrections to the 
effective action which carry derivatives. 

Our result means that SQED can be renormalized so that the inflationary 
production of scalars is not affected at one loop order. That is all we can 
conclude from the analysis of this paper, but the stochastic formalism can 
of course see to any order. It informs us that the quartic counterterm 5X 
can be chosen to cancel any (e^<y9*(y9)^ contribution to the effective potential 
H5] . This presumably implies that C2 = in (!99|) . However, there 
are unavoidable corrections to the effective potential at order {e'^(f*(f)^ and 
higher, which implies that the q's do not vanish for i > 3. So our result 
that Ci = is an artifact of being at one loop order which also happens 
to be repeated at two loop order. The nonperturbative outcome for SQED 
is that inflationary particle production engenders a nonzero photon mass, 
and the associated vacuum energy prevents the average scalar field strength 
from growing past ip*ip ~ iJ^/e^ [HI US]. The secular corrections to the scalar 
mode function (^Uh represent the scalar's initial response to the positive scalar 
mass associated with the vacuum energy of massive photons. 
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